Spin Dependent Collision of Ultracold Metastable Atoms 
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Spin-polarized metastable atoms of ultracold ytterbium are trapped at high density and their 
inelastic collisional properties are measured. We reveal that in collisions of Yb( 3 P 2 ) with Yb( 1 So) 
there is relatively weak inelastic loss, but with a significant spin-dependence consistent with Zeeman 
sublevel changes as being the dominant decay process. This is in strong contrast to our observations 
of Yb( 3 P 2 )-Yb( 3 P 2 ) collisional loss, which are, at low field, much more rapid and have essentially 
no spin dependence. Our results give a guideline to use the 3 P 2 states in many possible applications. 



Atoms with alkaline-earth-metal-like electronic struc- 
ture are under extensive study, partly due to their 
promise for a number of key applications. For example, 
the use of long-lived metastable P states (as well as the 
ground S state) has been explored as a useful quantum 
computing platform Also, the ultranarrow 1 S'o- 3 Po 

atomic resonance in a "magic wavelength" optical lattice 
may be highly competitive as a new optical frequency 
standard [6j . Finally, in the area of quantum simulation, 
there are several theoretical studies of the use of 3 Pj 
(J = 0, 2) atoms for studies of Hamiltonians with both 
spin and orbital degrees of freedom @, B , implementation 
of Abelian artificial g aug e potentials [9j , or simulation of 
Kondo lattice model |10j . 

Collisions of metastable P state atoms are interesting 
not only from a fundamental interactions viewpoint, and 
as the determining factor in collisional cooling schemes, 
but also as a crucial mechanism in several key applica- 
tions. Quantum gate phase imprinting via collisions has 
been proposed 0, [f| 



11] , as well as the exchange interac- 
tion in Kondo lattice model simulations flOj . Collisions 
can not only enable new physics but also inhibit desired 
applications. For example, in optical frequency standards 
collisional shifts can limit their accuracy, as well as lead 
to inelastic losses that destroy the atomic sample, pos- 
sibly on a time scale shorter than the desired interroga- 
tion time. Recent studies with metastable alkaline-earth- 
metal-like structure atoms showed inelastic collision rates 
as high as 10 -17 to 10 -16 m 3 /s in 3 Pj- 3 Pj collisions in 



Yb [12J and Sr [13j, |14| atoms. These rates, almost as 



high as the estimated elastic collisional rates, make the 
use of 3 Pj atoms difficult. 

The role of the collisions between the ground ( 1 5'o) 
and excited C*Pj) atoms is crucially important in many 
approaches 0, IE EH- In this respect, it is important 
to investigate the properties of collisions between the 
metastable triplet states and the ground state 14 , [lj| . 
In addition, the study of anisotropically interacting cold 
collisions is now a broadening area of stud y an d includes 
the rare-earth atoms and polar molecules ltjl- 24 1. 

In this Letter, we report the measurement of the in- 
clastic rate constants for Yb in both the 3 Po- 3 P2 and 



1 S'o- 3 P2 systems at several magnetic fields below 1 gauss. 
In particular, we observe strong spin dependence in the 
inelastic rates for Yb( 3 p2) m the 1 S'o- 3 p2 collisional sys- 
tem, with higher energy Zeeman sublevels having higher 
inelastic rates. This strongly suggests that the inelas- 
tic loss is dominated by Zeeman sublevel changing pro- 
cesses ( "ro-changing collisions" ) . These may be induced 
by the Landau-Zcncr transition between the entrance s- 
wave channel and the higher partial waves with lower 

f. In contrast, in the 3 p2- 3 P2 



magnetic sublevels [25 



26ij . in eoiiiraai, m me j~2~ i~2 
collisions we observe a much higher inelastic rate that is 
spin-independent. This is consistent with fine structure 
changing processes ( " J-changing collisions" ) or principal 
quantum number changing ("PQNC") processes. Our 
results represent a new, detailed study of 3 p2 physics, 
and also provide a roadmap for the use of the important 
3 P~2 state and challenge theory to provide a quantitative 
explanation of this anisotropic collision physics. 

In a previous experiment [12J, we studied 3 p2 collisions 
via indirect excitation of the 3 P>2 state to 3 P2. This led 
to a spin-unpolarizcd sample at the relatively high tem- 
perature of around 40 fiK. Here we prepare spin polarized 
samples at a much lower temperature, below 1 fiK, as de- 
scribed below. The experimental procedure is as follows. 
Yb atoms in a thermal beam generated from an oven at 
375°C are decelerated by a Zeeman slower with a strong 
1 So- 1 Pi transition at 399 nm, and then are loaded into 
the intercombination ( 1 5q- 3 Pi) MOT at 556 nm [27I. I28T] . 
FigurcQJb)-(d) show corresponding energy level diagram 
of Yb, timing diagram of the experiment, and schematic 
diagram of experimental setup, respectively. The atoms 
are transferred from the MOT to a crossed FORT at 
532 nm and evaporativcly cooled. The number of atoms 
is typically 3 x 10 5 . The BEC transition temperature 
(T c ) of our trap is 400 nK. All the work reported in this 
paper is done at 480 nK, just above T c , in order to main- 
tain a stable atom number. Our spin-selective excitation 
employs the very narrow linewidth 1 S'o- 3 P2 transition 
at 507 nm [29|, [3(j. The excitation efficiency to 3 P2 is 
about 10%, typically leaving about 2 x 10 4 3 p2 atoms in 
the trap. During the excitation wc apply a small bias 
magnetic field P>bias to spcctroscopically split magnetic 
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FIG. 1. (color online) (a) Schematics of the collision measure- 
ment. We studied both 3 P2- 3 Pi and 1 <S'o- 3 P2 collisional prop- 
erties, (b) Energy level diagram of ytterbium, (c) Timing di- 
agram for measurement of collision dynamics of 3 p2- 3 P2 and 
1 So~ 3 -P2 collision system. For measurement of 3 P2- 3 P2 col- 
lision system, unexcited 1 So atoms are immediately removed 
from the trap by applying a strong 399 nm laser. For mea- 
surement of 1 Srj- 3 P2 collision, the blast laser is irradiated af- 
ter interaction time, (d) Schematic diagram of experimental 
setup. Directions of 399 blast and probe lasers (not shown) 
are same as that of 507 nm excitation laser. 



sublevels mj of 3 P2 state. The applied magnetic field 
varies from 215 mG to 848 mG, corresponding to the 
Zeeman splitting of 0.45 MHz to 1.8 MHz. The excita- 
tion laser linewidth is measured to be less than 1 kHz and 
Dopplcr broadening is about 20 kHz. This is much less 
than the Zeeman splitting of 450 kHz at our lowest mag- 
netic field of 215 mG, thus ensuring selective creation of 
spin-polarized samples of atoms in a single Zeeman sub- 
level. 

We perform study of two different collisional systems, 
3 P 2 - 3 P 2 and 1 S - 3 P 2 . For the former, the leftover ^0 
atoms are immediately removed by strong excitation at 
399 nm (i.e. blast laser) [Fig. [He)]. After interaction 
time tint (where the 3 P 2 atoms collide with other atoms 
and undergoes inelastic loss), 3 P 2 atoms are repumped 
back to the 1 Sq state by 770 nm and 649 nm rcpumping 
lasers which are resonant to the 3 P 2 - 3 Si and 3 Po~ 3 Si 
transitions. A few milliseconds is required for complete 
rcpumping. Finally the number of repumped (^Sq) atoms 
are measured by absorption imaging using the 1 Sq- 1 Pi 
transition. For the latter (i.e. 1 Sq- 3 P 2 collision measure- 
ment), we simply leave the leftover 1 Sq atoms in the trap. 
Then the blast laser is irradiated after interaction time: 
just before the repumping. 

Inelastic atom-atom collisions in various channels are 
clearly observed for all Zeeman states. A selection of 
these state-dependent decay curves (those at highest and 
lowest magnetic field) are presented in Fig. [2J Figure 



[2^ (b) shows the time evolution of the number of 3 P 2 
atoms without (with) the So atoms. Thus the decay of 
atoms in Fig. is due only to 3 P 2 - 3 P 2 collisions. The 
measurements reveal a very high inelastic loss rate (ap- 
proximately the same as the estimated elastic rate) that 
is essentially independent of both mj and magnetic field 
strength. This behavior is consistent with the scenario 
that the dominant decay process is J-changing, which 



was inferred by our previous measurement [12(. It may 
also be due to PQNC collisions, as described in Ref. [131 ] . 

The decay of 3 P 2 state atoms in collisions dominantly 
with 1 S'o atoms is shown in Fig. [2Jd. It is important to 
note that in these decay curves the number of 1 Sq atoms 
is ten times the number of 3 P 2 atoms in either Fig. [2^ 
or [2Jd. Thus, although the decay curves are steeper in 
this raw data, as will be explained in detail later and 
shown in Fig. [3l 1 So- 3 P 2 collisions are less inelastic than 
3 P 2 - 3 P 2 collisions. One can easily recognize in Fig. [2Jd 
that the 1 So- 3 P 2 collisional decay behavior is markedly 
different from that of 3 P 2 - 3 P 2 collisions. The observed 
decay curves are strongly state-dependent. In particu- 
lar, the atoms in higher Zeeman levels in the 3 P 2 state 
(toj > 0) show a stronger field dependence in the decay. 
This clearly suggests the important role of m-changing 
collisions in the 1 S'o- 3 P2 collision. 

We extract 2-body inelastic loss rate coefficients from 
the observed decay rates using the following method. The 
decay is modeled by coupled differential equations: 



Tl e f3g e 7lg7l e f3 ee 7l e 
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(1) 



where n denotes the local density of atoms, subscript g 
(e) denotes ground (excited) state, (3 ge (/3 ee ) is the two 
body loss coefficient for ground-excited (excited-excited) 
state collisions. The trap lifetime for the ground state 
atoms is measured to be 30 s, thus we neglect back- 
ground collision loss because it is much longer than the 
present two-body collision loss. For the 3 P 2 - 3 P 2 colli- 
sion, Eqs. ([T]) become simple and are solved analytically, 
since all atoms in the ground state are removed from the 
trap(n g = 0). By spatially integrating the equation, we 
obtain the evolution in atom number 
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Here No e is the initial excited state atom number and 
Gee = (V 2e /Vi e )f3 e e- The effective volume is defined by 



V a . 



d 3 r[n e {r)/n 0c 



(3) 



where q is an integer number, nn e is the peak density 
of excited state atoms, and n e (r) is the spatial density 
distribution. 

The spatial density distribution depends on the dimen- 
sionlcss parameter 77: the ratio of the trap depth e t to 
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(a) 3 P 2 - 3 P 2 collision system 
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(b) 1 S - 3 P 2 collision dominant system 
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FIG. 2. (color online) Observed decay of trapped 3 p2 atoms as a function of time, (a) Time evolution of the number of 3 p2 
atoms without 1 So atoms. Decay rates of different magnetic sublevels mj at different magnetic field strengths are essentially 
the same, (b) Time evolution of the number of 3 P2 atoms with So atoms. In this case the number of 1 So atoms are 10 times 
larger than that of 3 P2 atoms. Thus the 3 P2 atoms dominantly collides with 1 So atoms. The behavior is markedly different 
from that of 3 P2- 3 P2 collision system. Decay rates are strongly spin-dependent. Note that we measured decay for all magnetic 
components in each bias magnetic field of 215, 307, 407, 596, and 848 mG. Part of data is shown in these graphs. 



the sample temperature fcgT in energy unit. In our ex- 
periment, 77 is 10 for the ground state atoms. The trap 
depth for each mj component in 3 p2 state strongly de- 
pends on the direction of external magnetic field and 
the polarization of FORT laser [30|, 31 1. To determine 
77 of the 3 P2 state, we measure the Stark shift for all 
mj states as a function of horizontal FORT laser power 
in various bias magnetic fields studied here. The ra- 
tio U e /Ug, where t/ e ( g ) is the trap potential of the 3 P2 
( 1 S'o) state, varies from 1.02 to 1.31. From the mea- 
surement, 77 for all mj states are calculated to be 7.2 
to 11.7. Since 77 keeps more than 7 for any mj com- 
ponents, large- 77 approximation can be used to calculate 
the effective volume. In the large-77 limit, n(r) is well 
approximated by the thermal density distribution [32j ]. 
n e (r) = n 0e exp[— C/(r)/fcsT], where U(r) is the trap po- 
tential. By approximating FORT potential to a trun- 
cated harmonic trap U(r) = e t (r / Rq) 2 Q(Ro — r), the 
effective volume can be written as 

3/2 



V q = i? 5 



<Pl 



(4) 



where 0(x) is the Heaviside step function, U(Rq) = et is 
the trap depth, and Rq is the boundary. 

We can extract the (3 ge coefficient by analyzing the de- 
cay of 3 P2 atoms in the presence of 1 So atoms, shown in 
Fig.[U(b). The data involves the 1 S U - 3 P2 collision as well 
as the 3 P2- 3 P2- Thus, the analysis is not so simple as Eq. 
©, because Eqs. ((T|) cannot in general be solved ana- 
lytically. However, if we assume depletion of the number 
of ground state atoms is negligibly small, Eqs. dJ) can 
be solved analytically. In this case, we obtain the atom 
number evolution in a large-77 approximation 

exp(-G ge N 0g t) 
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where G ge = f3 ge / [V{ 
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and N 0g is the ini- 



tial ground state atom number. Note that we can repro- 
duce Eq. © from Eq. {5J if we set N 0g = 0. 

We first fit the 3 P 2 - 3 P 2 collision data by using Eq. © 
and calculate /3™. Then decay data of 1 5'o- 3 P2 is fitted 
by using Eq. (|S|) with fixed As mentioned in Ref. 
[331 ] , the observed two-body decay rate j3 includes inelas- 
tic collision loss /3 ln and evaporation /3 cl due to clastic 
collision; i.e., j3 = (3 m + f/3 cl , where / represents the 
fraction of elastic collisions. The inelastic collision rate 
can be expressed as /3 m = (3/{fj +1), where 7 is the ra- 
tio of elastic to inelastic cross sections, 7 = a cl /a ul [33j . 
In a large-77 approximation, / can be approximated to a 
simple analytic equation [34| . In the present condition, 
/3 in //3 varies from 0.87 to 0.93. 

The inelastic collision rates for all mj states in various 
external magnetic fields are plotted in Fig. [3l Essentially 
no spin-dependence in the collision of two metastable 3 p2 
states is apparent, which is represented by triangles in 
Fig. |3l In particular, the spin states of mj > —2 show 
almost the same inelastic collision rates with the lowest 
energy state of mj — —2, which should not suffer from 
the m-changing collision. This behavior is quite different 
from the significant spin dependence theoretically pre- 
dicted in the collision of two metastable 3 p2 state at a 
high magnetic field 25 , 2f| . Therefore, our result clearly 
shows that the m-changing collision is not observed for 
all spin states at this low magnetic field and low tem- 
perature. Since the most likely decay process for the 
mj = —2 state is the J-changing collision from the 3 P2 
state to 3 Pn or 3 Pi, this would be dominant. PQNC col- 
lision may also contribute to the overall decay. However 
it is difficult to distinguish J-changing from PQNC col- 
lision unless we measure (3 ee for 3 Prj- 3 Pn collision, a task 
beyond the scope of this Letter. 

Note that the average of the obtained inelastic collision 
rate coefficients (3 ee of 4 x 10~ 17 m 3 /s is a little higher 
than the value of TOx 10 -17 m 3 /s, which is previously ob- 
tained at a 100 times higher temperature 12(. This inter- 
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csting temperature dependence is open to a further the- 
oretical investigation. In contrast, the significant spin- 
and field-dependence in the collision of the metastable 
3 P 2 state with the ground state 1 S is observed, which 
is represented by squares in Fig. [3] Especially, the spin 
states with higher energies (mj > —2) show higher in- 
elastic collision rates compared with the lowest energy 
state of mj = —2. Since the mj = —2 state differs from 
other spin states only in that the mj = —2 state does 
not suffer from the m-changing collision, it is natural to 
think that the dominant decay process is a m-changing 
collision in this 3 P2~ 1 <5'o collision at a low field. The quan- 
titative theoretical explanation of the detailed behaviors 
is an interesting future work [35| . 

In addition, we can claim that the mj = —2 state 
is rather stable against the collision with the 1 S'o atom. 
Since the decay of the 3 P 2 ijn.j = —2) atoms is dom- 
inated by the 3 P 2 - 3 P 2 collision even in the presence of 
the 1 5o atoms, as is shown in Fig. [3l it is quite difficult to 
accurately extract the /3 ge coefficient for mj = — 2 state. 
Although our analysis results in the value of /3 ge on the 
order of 10 -19 m 3 /s for the mj = —2 state, it may be 
much lower. In fact, our recent measurement on 1 Sq- 3 P2 
(m.j = —2) atoms in a tightly confined three dimensional 
optical lattice indicates that the upper limit of (3 ge is on 
the order of 10~ 20 m 3 /s [36| . This is promising to use in 
many applications. 

In conclusion, we have experimentally investigated col- 
lisional properties of spin polarized metastable 3 P 2 states 
of Yb Atoms. We reveal the significant spin-dependence 
in the collision of the metastable 3 P 2 state with the 
ground state 1 Sq, which strongly suggests that the dom- 
inant decay process is a m-changing collision. On the 
contrary, we observe essentially no spin-dependence in 
the collision of two metastable 3 P 2 states, which is con- 
sistent with the J-changing collision as a dominant decay 
process. Our results will trigger theoretical efforts to clar- 
ify these behaviors quantitatively and give a guideline to 
use the 3 P 2 states in many possible applications. 
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FIG. 3. (color online) Inelastic decay rate for all mj states in various external magnetic fields. Red triangles show inelastic 
decay rate /3™ of 3 P2- 3 P2 collision system. Essentially no spin-dependence is observed. Black squares show those of 1 So~ 3 P2 
collision system (/?!,")■ Significant spin- and magnetic-field-dependence is observed. 
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